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Abstract 

Phosphoryl transfer reactions figure prominently in energy metabolism, signaling, transport and motility. Prior detailed 
studies of selected systems have highlighted mechanistic features that distinguish different phosphoryl transfer enzymes. 
Here, a top-down approach is developed for comparing statistically the active site configurations between populations of 
diverse structures in the Protein Data Bank, and it reveals patterns of hydrogen bonding that transcend enzyme families. 
Through analysis of large samples of structures, insights are drawn at a level of detail exceeding the experimental precision 
of an individual structure. In phosphagen kinases, for example, hydrogen bonds with the 0 3p of the nucleotide substrate are 
revealed as analogous to those in unrelated G proteins. In G proteins and other enzymes, interactions with 0 3 p have been 
understood in terms of electrostatic favoring of the transition state. Ground state quantum mechanical calculations on 
model compounds show that the active site interactions highlighted in our database analysis can affect substrate 
phosphate charge and bond length, in ways that are consistent with prior experimental observations, by modulating 
hyperconjugative orbital interactions that weaken the scissile bond. Testing experimentally the inference about the 
importance of 0 3p interactions in phosphagen kinases, mutation of arginine kinase Arg 28 o decreases k cat , as predicted, with 
little impact upon K M . 
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Introduction 

Enzymes that catalyze the transfer of a phosphate from ATP are 
widespread in biology. Free energy liberated in hydrolysis of 
phosphoandyride bonds, such as those in ATP, is essential for 
cellular energy metabolism, motility and the generation of 
transmembrane potentials. Much has been learned about 
enzyme-catalyzed phosphoryl transfer through the detailed study 
of individual enzymes. The most frequently cited mechanisms 
include; a) precise positioning of substrates for phosphate transfer; 
b) base-assisted activation of the nucleophile; c) activation of the 
electrophile; and finally d) electrostatic stabilization of the 
transition state [1]. Additional mechanisms have been proposed 
for specific enzymes, including: proton relay systems [2], substrate- 
assisted catalysis [3] and strain of the P and y phosphate groups 
[1]. Phosphoryl transfer enzymes are a broad class within which 
several mechanisms may have evolved. However, one is struck by 
the diversity in mechanistic proposals and the lack of consensus on 
key characteristics of active sites that might implicate common 
elements of mechanism that might bridge across diverse enzyme 
families. 



The Protein Data Bank (PDB) provides an opportunity for a 
top-down analysis of common active site configurations, including 
hydrogen bonding interactions with water. In this work, 
phosphoryl transferases are used to test the postulate that 
comparative analysis can complement traditional reductionist 
investigations of individual systems, providing additional enzymo- 
logical insights. Three classes of structures are compared here: 
those that cleave an — P y bond (Figure la); those that cleave a 
P a — 0 3a bond (Figure lb), and a "control" group of those that 
bind nucleotide without catalyzing phosphoryl transfer (Figure lc). 
This work stems from the premise that important atomic 
interactions may have been hidden by the experimental error in 
the analysis of individual structures. We explore whether they can 
be revealed through statistical analysis of populations of structures 
of proteins sharing similar function, examining the spatial 
distributions of atoms at a local level. Our statistical analyses 
reveal that interactions with the bridging oxygen of a scissile 
phosphoanhydride bond, previously noted for select enzymes 
[2,4,5,6], are present in diverse enzymes. Interactions with non- 
bridging (3-phosphate oxygens are also widespread in active sites. 
The analysis is put to brief experimental test through kinetic 
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analysis of an active site mutant affecting interactions with the 
bridging (3-oxygen in the reaction of arginine kinase, an enzyme in 
which such interactions had not previously been implicated. 
Finally, quantum mechanical calculations reveal a possible 
stereochemical rationale for the observed interactions with both 
bridging and non-bridging (3-oxygens, adding to a number of 
plausible proposals that seek to explain the bond selectivity and 
mechanism of phosphoryl transfer. 

Materials and Methods 

Database analysis 

Structure sets. Coordinates for protein-nucleotide complex- 
es (Table SI) were downloaded from the PDB, including those 
containing ATP, UTP, TTP, GTP, CTP or their analogs. Two 
sets of data were compiled, one containing structures up to 2.0 A 
resolution, where solvent molecules were defined with confidence, 
and a larger group (inclusive of the first) up to 2.7 A resolution 
from which solvent water was excluded, due to insufficient 
confidence in their positions at this resolution. Of the 1866 
protein-NTP complexes to 2.7 A, as of October 3 rd , 2012, 305 
were non-redundant and satisfied inclusion criteria (Table S2) 
that, for example, excluded structures in conformations known to 
be non-catalytic. Of the 305 vetted structures, 134 were at 
resolutions higher than 2.0 A. Molprobity's Reduce [7] was used 
to flip asparagine, glutamine and histidine side chains for optimal 
hydrogen bonding. 

Structures were grouped by the site of bond cleavage in the 
ligand: 1) at the 0 3 p — P y bond (Figure la); 2) at the P a — O 3oc 
(Figure lb); and 3) those in which the ligand binding site is non- 
catalytic (Figure lc). The groups had 155, 100, and 50 PDB 
structures in the 2.7 A set, respectively and 72, 45 and 17 
structures in the 2.0 A set. An underlying premise was that 
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Figure 1. Enzyme structures can be categorized according to 
the fate of the bound nucleoside triphosphate (NTP). A) 

Phosphoryl transfer in which the 0 3 p_P Y bond is cleaved. B) Reactions 
in which the Po_0 3a bond is cleaved and C) Structures where the bound 
NTP does not undergo a chemical reaction. Red lettering indicates the 
atoms in the scissile bond and red arrows depict the transfer of 
electrons in going from reactants to products. 
doi:1 0.1 371 /journal.pone.01 0831 0.g001 



enzymes might share catalytic features regardless of base-type 
(GTP, ATP, etc.). The set of proteins with non-catalytic binding 
sites serves as a control to distinguish potentially catalytic 
interactions from those primarily involved in ligand-binding. 
Comparisons between the sets of enzymes with P a — C>3 a and 
0 3 p — P Y cleavage sites also help to distinguish binding and 
catalytic interactions, due to spatial separation of binding and 
catalytic sub-sites. 

Identification of interactions. Coordinates were expanded 
according to the crystallographic symmetry and riding hydrogen 
atoms were added. These hydrogens had not been seen directly by 
x-ray crystallography, so their positions were inferred from the 
heavy atom coordinates using Reduce from the MolProbity 
program [7]. Assessed interactions included hydrogen bonds and 
coordination of metal ions, all analyzed using an in-house Python 
program. 

Heavy atom distance and angle criteria were based on previous 
work [8]: O—D distance <3.5 A, P— O—D angle >90°, where D 
is the hydrogen bond donor, P is the ligand phosphorus atom, and 
O is the oxygen atom of the ligand. An additional criterion was 
added: 0"D — R angle ^90° where R is the antecedent atom of a 
hydrogen bond donor. These criteria depend only on the heavy 
atoms directly observable and can therefore be applied directly to 
medium resolution crystal structures. Assessments were also made 
using additional criteria possible once riding hydrogens were 
added: 0—H<2.7 A, P— O—H>80° and O—H— D>80°. In- 
spection revealed that a small number of atom pairs satisfied the 
criteria without interacting directly, because they were brought 
into proximity by coordinating with a common metal ion. These 
were excluded using criteria shown in Figure SI. Metal 
interactions were subject only to a distance cut-off of 2.8 A 
derived from analysis of metal coordination in the Cambridge 
Structural Database. 

Classification and comparison of interactions. Ligand 
oxygens were organized into six groups: Non-bridging y-, (3- and 
cx-oxygens, and bridging (3-, oc- and 05 '-oxygens. Distributions of 
enzyme-ligand interactions for each of the groups were compared 
using pair-wise two-tailed t-tests in RStudio using pooled standard 
deviations. When comparing the numbers of interactions with 
different atoms for the different structure sets, the Benjamini & 
Hochberg correction for multiple comparisons was applied. For 
each comparison, the null hypothesis was tested that the two 
sample means could be drawn from the same population, at a 
significance level of oc = 0.05. Results are fully tabulated in the 
supplemental material (Table S3 and Table S4). 

Quantum mechanical analysis of interactions 

The Gaussian 09 program [9] was used for all electronic 
structures calculated with the B3LYP/6-31 1++G(d,p) level of 
theory [10]. Bulk solvent was modeled implicitly with the 
polarizable continuum model (PCM) using a solvent dielectric of 
78.36 for water [11]. This level of theory had previously been 
shown to be appropriate for capturing solvation effects [12]. 
Natural bond orbital (NBO) analysis [13] was performed with the 
GenNBO5.0 program [14] using the B3LYP/6-31 1++G(d,p)- 
optimized structures and HF/cc-PVTZ model chemistry [15]. 
The NBO program facilitates a Lewis-like description of electron 
density by transforming nonorthogonal atomic orbitals from a HF 
wave function into natural bond orbitals. Second order perturba- 
tion E(2) energy, given by Equation 1, estimates the magnitude of 
electron derealization and hyperconjugation. 
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E(2) = q t 



Ej-Et 



qi is the donor orbital occupancy, and Ej - Ei is the energy gap 
between donor orbital i and acceptor orbital, j. Fy is the Fock 
matrix element which describes donor and acceptor orbital 
overlap. In this work, the donor will be a lone pair orbital of the 
y-phosphate oxygens and the acceptor orbital is G*(0 3 |3 — P y ), the 
antibonding orbital of the 0 3 p — P y bond in methyl triphosphate. 
Methyl triphosphate has been used previously, to model interac- 
tions in ATP (Figure 2) [16]. 

The influence of a N — H« M 0 3 (3 hydrogen bond on the 0 3 p — P y 
bond length and hyperconjugation was evaluated with Structures 
1 through 6 (Figure 2a— f). (Similar results were obtained with the 
addition of Mg ++ coordinated to (3- and y-oxygens.) Structure 7 
was used to evaluate the effects of a hydrogen bond to a 
nonbridging (3-oxygen (Figure 2g). N-methylacetamide was used 
to model a peptide backbone (neutral) hydrogen bond donor. The 
N — H«"Op hydrogen bond was weakened by increasing its length 
from 2.5 A to 3.5 A in 0.1 A increments for Structures 1 through 
3 and in 0.2 A increments for Structures 4 through 7. These 
changes were made by translating the N-methylacetamide, fixing 
the 0 3 p and N atoms at a specific hydrogen bond distance, then 
optimizing other aspects of the structure with the N — H»»»Op 
angle fixed at 170°. Structures 2 through 6 were used to probe the 



electronic effects of either protontion or hydrogen bonding at the 
y-oxygens of methyl triphosphate (Figure 2b— f). For protonated 
(1) structures the 0 3 p — P y — O y — H dihedral angle was fixed at 180°. 
Acetamide was used as a model for an asparagine side chain in 
Structure 5, and 1 -propylaminium was used as a model for a lysine 
side chain in Structure 6. All hydrogen bonds to y-oxygens, were 
fixed at a D«"O y distance of 2.8 A, where D is the hydrogen bond 
donor heavy atom. The D«"O y — P y — 0 3 p and the (D)H«"O y — 
P y — 0 3 p dihedral angles were each fixed at 180°. Changes in 
0 3 p — P y bond lengths and in E(2) energies were calculated relative 
to each molecule without a hydrogen bond between N- 
methylacetamide and Op. 

Mutagenesis and Kinetics 

Mutants of Limulus polyphemus arginine kinase (AK) were 
made using the QuikChange mutagenesis kit (Stratagene, Inc.). 
Expression, purification and kinetic analysis of mutants R280L 
and R280K followed previously published methods [17], though 
only R280K had sufficient activity for full analysis. Initial velocities 
were measured using a 6x6 grid of substrate concentrations: ATP 
at 0.2, 0.4, 0.8, 1.2, 1.6, and 2.0 mM, and arginine at 0.2, 0.4, 0.8, 
1.6, 3.2, and 7.0 mM. Duplicate measurements were made with 
separate protein preparations. Steady state parameters K ia , K m 
and V max for a random order sequential bimolecular-bimolecular 
(bi-bi) reaction mechanism were fitted to the data by nonlinear 
least squares using SigmaPlot. 
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Figure 2. Models used to test the dependence of hyperconjugation and 0 3p _P y bond length on enzyme-ligand interactions. N- 

methylacetamide (A) was used to model a (neutral) protein backbone amide hydrogen bond to 0 3(3 , using methyl triphosphate to model an NTP 
nucleotide. Structures 2 (B) and 3 (C) were used to investigate the effects of protonation at the y-oxygens. Additional active site hydrogen bonds, 
represented in Structures 4 through 6 (D-F), were used to assess the secondary effects of different types of 0 Y hydrogen bonds. Acetamide (E) 1- 
propylaminium (F) were used to model asparagine and lysine side chains respectively. Structure 7 (G) was used to investigate the impact of hydrogen 
bonding at a nonbridging p-oxygen on hyperconjugation and 0 3 p_P Y bond length. Hydrogen bonds are shown with dashed red lines. 
White = hydrogen, gray = carbon, red = oxygen, blue = nitrogen, orange = phosphorus. 
doi:1 0.1 371 /journal. pone.01 0831 0.g002 
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Results & Discussion 

Database analysis of enzyme-nucleotide interactions 

Nucleotide interactions included hydrogen bonds with the 
enzyme, water or other substrates, as well as interactions with 
metal ions. The numbers of interactions were compared between 
phosphoryl transfer enzymes (Figure la) and two comparison 
groups (Figure lb, c) to distinguish potentially catalytic from non- 
catalytic associations. T-tests reveal two overall trends (Table S3 
and Table S4), that interactions with both non-bridging (3-oxygens 
and with the bridging oxygen of a scissile phosphoanhydride bond 
are more prevalent in phosphoanhydride bond-cleaving sites than 
in non-catalytic nucleotide-binding sites. 

For non-bridging (3-oxygens, the mean numbers of interactions 
in both 0 3 p — P y -cleaving sites (|i = 4.2, a = 0.2) and P a — O 3oc - 
cleaving sites (|i=3.3, a = 0.1) are significantly greater than in 
non-catalytic sites ((1 = 2.5, a = 0.2; p<0.05; 2.7 A; Figure 3a, 
Table S3a). The database analysis extends the finding in Ras, by 
Fourier Transform Infrared (FTIR) spectroscopy, of stronger 
interactions with (3-oxygens than (X- or y-oxgyens in [18], 
indicating that strong Op interactions are widespread among 
enzymes. Enzymes that cleave the 0 3 p — P y bond favor neutral H- 
bond partners with the non-bridging (3-oxygens, whereas positively 
charged hydrogen bond donors are common in non-catalytic NTP 
binding sites. Backbone amide Op interactions have been 
implicated with Walker A motif P-loops [1,19], so the pre- 
eminence of neutral hydrogen bonding in the database indicates 



that analogous interactions are favored in a more general array of 
phosphoryl transfer enzymes. 

Although not every protein interacts with the bridging 0 3 p, 
interactions are about twice as common in 0 3 p — P Y cleaving sites 
(mean number of 0.8, a = 0.1) as in either P a — O 3oc cleaving sites 
(|i = 0.4, a = 0.1) or non-catalytic NTP-binding sites ((1 = 0.5, 
a = 0.1; p<0.05; 2.7 A; Figure 3b, Table S3a). (The statistics 
quoted are for the larger data set, but similar results are also 
obtained with the smaller 2.0 A set that included solvent waters, 
Table S3b and Table S4b). Hydrogen bonds with the 0 3 p most 
commonly involve neutral donors, particularly from backbone 
nitrogens (Table S5). Other donors include water and positively 
charged amino acids, especially arginine. 

Interactions with bridging 0 3 p have previously been cited for G 
proteins and several kinases, including 6-phosphofructo-2-kinase/ 
fructose-2,6-fosphosphatase and pyruvate kinase [2,5,20]. Thus, 
the structural survey captures results expected from detailed prior 
analyses and shows that analogous interactions actually occur in a 
diverse array of phosphoryl transfer enzymes including the MutL 
DNA mismatch repair protein [21], the DEAD-box RNA helicase 
[22], and Phosphoenolpyruvate carboxykinase [23]. The donor is 
often, but need not be, within some type of conserved motif, for 
example at position 4 of the Walker A motif (p-loop; Table S5). 
0 3 p interactions have previously been rationalized in terms of 
stabilizing development of negative charge on the leaving group in 
proposed dissociative transfer mechanisms [2,5,20], though 
Warshel has noted that associative mechanisms would likely also 
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Figure 3. Mean number of enzyme-ligand interactions. Interactions are shown with nonbridging (3- (A), bridging (3- (B) and bridging a-oxygens 
(C) in the 0 3) 3_P Y cleaving-, Po_0 3a cleaving- and non-catalytic-NTP-binding sites. The resolution cutoff of structures used is 2.7 A, and water is not 
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benefit from electrostatic stabilization of increased electron density 
on the [24,25]. Our survey demonstrates O3P interactions 

occur in diverse enzyme families, such as the phosphagen kinases, 
N-Acetyl-L-Glutamate kinase and Adenosine kinase, in which the 
mechanisms have been reported to be significantly associative 
[26,27,28,29]. Thus, the interaction may transcend the associa- 
tive/dissociative mechanistic categorization. The expanded set of 
enzymes exhibiting 0 3 p interactions represent a wide array of 
protein folds and differ in other aspects of active site configuration. 
Indeed, the proposed primary catalytic mechanisms of these 
enzymes are quite diverse, so it was not obvious that such a broad 
array of enzymes would share analogous interactions. The much 
greater incidence of 0 3 p interactions among enzymes cleaving at 
0 3 p — P y over P a — C>3 a , and over proteins with non-catalytic 
nucleotide binding sites, suggests that most of these interactions 
have evolved to have a role in the chemistry and not just binding. 

Interactions with C>3 a in enzymes that cleave the P a — 03 a bond 
are more frequent than in the comparison groups, but more 
marginally than the 0 3 p interactions. Only the difference between 
P a — 03 a cleaving sites ((1 = 0.4, a = 0.1) and O3P — P y cleaving 
sites (|i = 0.3, G<0.1) reaches statistical significance (p<0.05, 
Table S3 a). 

A preference was found for positively charged hydrogen bond 
donors in non-catalytic binding sites relative to the catalytic sites 
(0 3 p — P Y and the P a — 0 3a combined group; Figure 4a, b). 
Differences manifested particularly at the non-bridging y-oxygens 
(catalytic: (1 = 2.0, a = 0.2; non-catalytic: |i= 1.4, a = 0.1; p<0.05; 
2.7 A; Table S6)). The number of neutral hydrogen bond donors 
in catalytic sites exceeded the average in non-catalytic sites 
(Figure 4c, d) particularly for the non-bridging (catalytic: (1 = 0.8, 
a = 0.2; non-catalytic: (1= 1.9, a = 0.1) and bridging (3- oxygens 
(catalytic: 0.2, a = 0.1; non-catalytic: (1 = 0.4, a<0.1; p<0.05; 
2.7 A; Table S6). There are several possible explanations for the 
observed differences in disposition of charged interactions, 
including the possibility that non-catalytic binding sites have 
evolved without pressure for efficient dissociation of reaction 
products. 

Quantum mechanical modeling of nucleotide 
interactions 

Several mechanistic rationales have been offered for the 0 3 p 
interactions in G proteins and related enzymes 
[2,4,5,6,20,25,28,30,31,32]. This section examines if such inter- 
actions could (also) affect nucleotide stereoelectronics in light of its 
recently appreciated importance in the modulation of phosphoan- 
hydride bond energies [12,16]. The total of hyperconjugative 
effects has been estimated at —300 kj/mol, with hydrogen bonds 
to y-oxygens each lessening this by ~20 kj/mol [16], so 
modulation of active site interactions could plausibly impact 
reactions with activation barriers like the 34 kj/ mol of arginine 
kinase [17]. Our approach is Natural Bond Order analysis [13] 
which can inform us of the ground state disposition, but is not 
applicable to the transition structure, so the approach allows only 
limited inferences about rate enhancement. 

For G proteins, the O3P interactions have been rationalized 
through electrostatics [5]. Hydrogen bonding or cation interac- 
tions would favor charge accumulation on O3P in the progression 
towards a dissociative transition state, or it could favor the 
transient accumulation of electron density in an associative 
mechanism [24,25]. It has also been noted that favoring charge 
accumulation on 0 3 p would also tend to increase the susceptibility 
of P y (and Pp) to nucleophilic attack [28]. A dissociative 
mechanism has been cited as consistent with FTIR experiments 
on Ras-GAP [33] that indicated reduction of (3- and increase in y- 



phosphoryl group bond orders. It is thought that electrostatic 
stabilization of negative charge at the (3-oxygens draws negative 
charge (by induction) out of the 0 3 p — P y bond onto the 0 3 p 
oxygen, which is then partially transferred to the non-bridging (3- 
oxygens through resonance [33,34,35]. A small molecule quantum 
study of a series of organic phosphates indicated, however, that 
hyperconjugation was likely to be a greater factor than induction 
in scissile bond lengths [12]. Hyperconjugation is a stereoelec- 
tronic effect, where electron density is transferred from an electron 
rich donor orbital to the antibonding orbital of a neighboring 
bond [36]. An example is depicted in Figure 5, where electron 
density is transferred from a y-oxygen lone pair orbital in methyl 
triphosphate into the antibonding orbital of the 0 3 p — P y bond. 
Rationalization of the mechanistic impact of 0 3 p interactions has 
already been controversial [5,24,25,28]. The goal here was not to 
arbiter this discussion, but to examine the plausibility of 0 3 p 
interactions having impact mediated through hyperconjugation, 
possibly in addition to a mixture of other effects. 

Past quantum mechanical modeling of methyl triphosphate 
[12,16] has shown that enhanced hyperconjugation renders y- 
oxygens more positive and the 0 3 p more negative as electron 
density is transferred from the y-phosphoryl group to the G* 
0 3 p — Py antibonding orbital, lengthening the scissile bond. 
Qualitatively, the changes in charge and bond length are the 
same as those proposed for electrostatic 0 3 p interactions [5] and 
observed through FTIR on RAS-Gap [33]. Thus the task is not to 
rule out the possibility of electrostatic contributions to the 
mechanism but rather to determine whether hyperconjugation 
could be a major contributor as well. Whether interactions with 
the (3-oxygens can enhance hyperconjugation in a nucleotide, is 
tested through quantum mechanical calculations on a methyl 
triphosphate model system with a hydrogen bond to one of its (3- 
oxygens from N-methylacetamide (Figure 2). The dependence of 
hyperconjugation upon the N — H»»»03p hydrogen bond is 
investigated by calculation of the 0 3 p — P y bond length as the 
length of the N — H"«0 3 p hydrogen bond is perturbed. 

Density functional calculations on Structure 1 (Figure 2a) show 
that a single hydrogen bond, between N-methylacetamide and the 
O3P oxygen in methyl triphosphate, induces an increase in 
hyperconjugation and a —0.05 A elongation of the 0 3 p — P Y bond 
(Figure 6). This affirms the hypothesis that the 0 3 p — P y bond 
length can, at least in principle, be weakened by O3P hydrogen 
bonding acting through hyperconjugation. 

Two factors determine the magnitude of hyperconjugation; the 
degree of donor and acceptor orbital overlap, and the orbital 
energy gap. The orbital overlap (Fy in Eq. 1) and the orbital 
energies of n(O y ) and G*(03p — P y ) were calculated for Structure 1 
(Figure 2a), as a function of hydrogen bond length (Figure 7). 
There is a strong dependence of the energy of the G*(03p — P y ) 
antibonding orbital upon strength of the hydrogen bond 
(Figure 7 a). However, neither the energy of the non-bridging y- 
oxygen lone pair orbitals (Figure 7b) nor the orbital overlap 
(Figure 7c) are affected. As the hydrogen bond is shortened (made 
stronger), the orbital energy gap decreases, increasing hypercon- 
jugation and elongating the 0 3 p — P y bond. 

The relationship between 0 3 p — P y bond elongation and 
hyperconjugation was confirmed by protonating the y-oxygens. 
Protonation of the y-oxygens decreases the hyperconjugative 
interaction: n(O y )^G*(0 3 p — P y ) [12,16]. We now find that 
protonation also decreases the dependence of the 0 3 p — P Y bond 
length on the strength of the D"«0 3 p hydrogen bond (Figure 8), 
consistent with the impact of an O3P hydrogen bond upon O3P — 
Py elongation being significantly mediated by hyperconjugation. 
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Hyperconjugation is also subject to other "secondary" pertur- 
bations that might be of greater biochemical relevance; hydrogen 
bonds added at the y-oxygens interfere with hyperconjugation 
[16]. The primary dependence of — P y elongation on the 
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Figure 6. The 0 3 p — P y bond length and hyperconjugation 

increase with decreasing D«««0 3 p hydrogen bond length in 

Structured Calculated E(2) energy is of the n(0 Y )^a*(0 3 p— P Y ) 

hyperconjugative interaction. D is the N — H hydrogen bond donor in 

N-methylacetamide. 
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strength of the 0 3 p hydrogen bond was then recalculated, after the 
addition of (secondary) hydrogen bonds at y. Several y-oxygen 
hydrogen bond donors were considered: water, acetamide 
(mimicking an asparagine) and 1 -propylaminium (mimicking a 
lysine; Figure 2). Hydrogen bonds to y-oxygens were fixed at 
2.8 A, while the hydrogen bond between N-methylacetamide and 
the 0 3 |3 oxygen was varied from 2.5 A to 3.1 A in 0.2 A 
increments. 

Hydrogen bonds from neutral donors (water or acetamide) to a 
y-oxygen have little effect (Figure 9). Only the salt bridge bond 
from a positively charged 1 -propylaminium donor modulates the 
primary interactions. The effect is modest, reducing bond 
elongation by about 15%, but it is possible that multiple O y salt 
bridges could have an additive effect. 
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Figure 8. Impact of protonation of the y-oxygens. Protonation at 
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Quantum mechanical analysis of interactions with non-bridging 
oxygens - Of interest is whether hyperconjugation can explain the 
higher frequencies of interactions with non-bridging (3-oxygens in 
— P y cleaving active sites, relative to non-catalytic binding sites 
(Figure 3). Similar to our study of 0 3 p interactions, QM 
calculations were performed for N-methylacetamide hydrogen- 
bonded to a non-bridging (3-oxygen in Structure 7 (Figure 2g). 
Changes in the n(0 Y )^G*(0 3 p — P y ) hyperconjugative interaction 
were calculated as the hydrogen bond length was varied from 
2.5 A to 3.5 A in 0.2 A increments. 

The impact of a single hydrogen bond with a non-bridging (3- 
oxygen is smaller than that of an O3P hydrogen bond (Figure 10). 
However, the database analysis revealed an average of four 
enzyme interactions summed across the two non-bridging (3- 
oxygens (Figure 3a), nearly saturating the lone pairs. Multiple 
hydrogen bonds to a single oxygen in methyl triphosphate had an 
additive effect on hyperconjugation [16]. Thus the aggregate effect 
of interactions with non-bridging (3-oxygens could be commensu- 
rate with those of the 03(3. 
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In summary, hyperconj ligation weakens the scissile — P y 
bond. Hydrogen bonds to the phosphate oxygens can either 
amplify or reduce this effect, according to proximity with the 
scissile bond. The greatest weakening of the scissile bond comes 
from interactions with 0 3 p, proximal to the scissile bond, with 
lesser weakening from interactions at non-bridging Op that are 
two bonds removed. By contrast, hydrogen bonds to y-oxygens, 
vicinal to the scissile bond, interfere with the hyperconj ugation, 
strengthening the scissile bond [16]. Changes in hydrogen bonding 
are inherent to the act of substrate binding, but of interest is 
whether particular types of interactions have been selected in 
enzyme evolution. Here we test the hypothesis that enzymes 
enhance hyperconj ugation by maximizing binding to (3- and 
minimizing binding to y-oxygens. 

Noting that the desolvation of the non-bridging (3- and y- 
oxygens on substrate binding would have opposite effects upon 
0 3 (3_P y bond elongation, the 2 A structure database was queried 
for the fractional saturation of lone pairs with enzyme or solvent 
interactions, comparing the total five lone pairs of non-bridging (3- 
oxygens to the 8 of non-bridging y-oxygens. As predicted, the non- 
bridging (3-oxygens retain higher numbers of interactions than 
either the non-bridging a- or y-oxygens (Figure 1 1, Table S7), and 
there are more numerous interactions with the bridging than 
either the O 3oc or 0 5 '. The favoring of interactions with the (3- 
phosphoryl group that we see in the database analysis is consistent 
not only with the QM-based prediction, but with the FTIR studies 
of Ras that showed stronger binding of the (3-phosphoryl group 
[37]. 

Consistency of the observed disposition of active site interactions 
makes a circumstantial case for the importance to enzymes of 
hyperconjugation-mediated weakening of the 0 3 p — P y bond. The 
ground state QM calculations show clearly that bond selectivity is 
impacted, so enzymes could be taking advantage of hyperconjuga- 
tion in specifying which of the phosphoanhydride bonds is to be 
cleaved. Understanding whether modulation of hyperconjugation 
has a role in rate enhancement is a more challenging question 
about which we cannot be as definitive. Bypassing the question of 
rate enhancement, even a role in bond selectivity alone could be of 
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doi:1 0.1 371/journal.pone.01 0831 0.g01 0 

sufficient selective advantage to explain the prevalence of 
hyperconj ugation-mediating interactions in the structural data- 
base. 

Ground state quantum mechanical modeling provides, through 
hyperconjugation, a unifying rationalization for interactions with 
the phosphate oxygens. The calculated changes in bond order and 
partial charge are qualitatively consistent with those inferred from 
kinetic isotope and FTIR studies [33,34], and are consistent with 
the dispositions of (3- and y-oxygen interactions seen in the 
database analysis. Hydrogen bonding at the bridging 0 3 p and 
non-bridging (3-oxygens closes the energy orbital gap, increasing 
the transfer of electron density into the G* antibonding orbital and 
lengthening the 0 3 p_P y bond, while interactions at y-oxygens 
have the opposite effect. Thus, differences in the interactions with 
(3- and y-oxygens can modulate the strength of the scissile bond in 
a mechanism mediated through hyperconjugation. 

To affect rate, the activation barrier must be lowered by some 
combination of stabilizing the transition state (TS) or destabilizing 
the ground state. The (spontaneous) existence of hyperconjugation 
in isolated nucleotides means that it can be stabilizing the ground 
state. However, inter-molecular interactions that increase hyper- 
conjugation beyond its natural equilibrium level would destabilize 
the substrate. Yang and Cui examined whether the anomeric 
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effect in myosin destabilized the ground state [38]. A modest 
increase in ATP calculated potential energy, as the P y — 0 3 p bond 
was lengthened, was deemed insufficient for substantial rate 
enhancement. However, the argument implicitly depends on a 
reciprocal relationship between bond length and hyperconjuga- 
tion, such that distortions of the bond induce a change in 
hyperconjugation to the corresponding level of the enzyme-bound 
substrate, a questionable assumption. A definitive assessment of 
the influence of hyperconjugation on transition state stabilization, 
has been stymied by inapplicability of NBO calculations to the TS, 
and by flat energy surfaces (calculated for myosin and RasGAP 
reactions) from which it is difficult to distinguish different possible 
mechanisms and estimate entropic terms [25,38]. These energy 
calculations did, however, imply that the phosphoryl transfer 
reaction was mixed associative/dissociative in character. Natural 
resonance theory calculations [39] for model associative and 
dissociative reactions showed that the hyperconjugation increased 
on approach to a dissociative transition state, and remained near 
constant for associative reactions [40] . Thus, for the Ras-Gap and 
other enzymes with a mixed mechanism, overall interactions that 
enhance hyperconjugation are predicted to stabilize the TS. 
Therefore the modulations characterized in this work could 



plausibly be enhancing rate as well as bond selectivity. However, 
we lack a detailed understanding of how the interactions might 
change on progression to the TS, so caution is needed in 
speculating about rate enhancement. Furthermore, active site 
interactions commonly have multiple effects which are not 
comparatively evaluated in the current work. Thus, the current 
work adds hyperconjugation as one of the plausible rationaliza- 
tions of phosphoryl transfer rate enhancement, without ruling out 
the previously discussed electrostatic effects [5,24,25]. 

Experimental characterization of an arginine kinase 
mutant 

Among the diverse enzymes shown by the survey to have 0 3 p 
interactions, the phosphagen kinases are of particular interest. 
Physiologically, these enzymes buffer cellular ATP levels through a 
reversible transfer of the y-phosphoryl to a guanidine substrate 
such as creatine or arginine [41]. The enzymes have no structural 
homology to those where 0 3 p interactions were previously 
implicated, and have not been reported to share mechanisms of 
catalysis. Indeed, in contrast to the dissociative mechanisms 
postulated for G-proteins and other enzymes with previously noted 
0 3 p interactions [2,5,20], a more associative mechanism was 
suggested for phosphagen kinases [26,27]. A high (1.2 A) 
resolution crystal structure of an arginine kinase transition state 
analog complex added circumstantial evidence [42]. The nitrate, 
mimicking a transferring phosphoryl was situated between the 
guanidino N and nucleotide 0 3 p which were separated by 6 A, 
consistent with partial axial bonding of a pentavalent P Y , and the 
angles of approach for guanidino N and 0 3 p were within a few 
degrees of optimal for nucleophilic attack in an early stage of the 
reaction, before deprotonation of the guanidino N. This suggested 
an active site in which enzyme interactions were optimized for a 
significantly associative reaction. Such contrasts to the aforemen- 
tioned G proteins, together with the availability of high resolution 
structure, rendered arginine kinase particularly suitable as an 
experimental test of the predictive power of the database analysis. 
Below we report an experimental validation of the importance of 
0 3 p interactions highlighted by the database analysis, noting that 
the steady state kinetics can be used to investigate the overall role 
of an interaction, but not detailed mechanism. 

In the arginine kinase structure, the 0 3 p interacts with Arg 280 
and Arg 126 . Arg 12 6 interacts with both 0 3 p and O y , so Arg 280 , 
which also interacts with an cx-nonbridging oxygen of the 
nucleotide substrate (and an enzyme aspartate) was chosen for 
mutagenesis, because it did not have potentially confounding 
interactions with y-oxygens (Figure 12). Preliminary analysis 
showed that the activity of a R280K mutant was sufficient for 
steady-state kinetics analysis, but not that of a R280L mutant. 
Relative to the wild-type (WT), k cat of R280K was decreased from 
104 s 1 to Is \ whereas K M values for the two substrates 
remained within 3-fold of WT (R280K: K M (Arg) = 0.8 mM; 
K M (ATP) = 0.5 mM. WT K M (Arg) = 0.3 mM; 

K M (ATP) = 0.3 mM) [43]. The much larger impact upon rate 
indicates that the effect of the mutation is predominantly upon the 
chemical reaction as opposed to substrate binding. 

The arginine kinase R280K is analogous to the R291K 
mutant of the homologous enzyme, creatine kinase (GK) in 
rabbit muscle, which yielded similar kinetics (k cat = 2.2 s -1 and 
K M (ATP) = 1.6 mM) [44]. It was suggested that this mutation 
abrogated catalysis by misaligning the substrates. However, near 
wild-type K M values argue against gross distortions of substrate- 
binding, although the effects of subtle modulations of alignment 
cannot be ruled out [45,46]. Our database survey now highlights 
the analogy of the 0 3 p interaction with the corresponding ones in 
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Figure 12. Arginine kinase active site mutation. Selected 
interactions of the nucleotide in the crystal structure of the transition 
state analog complex of Horseshoe Crab Arginine Kinase (AK, PDB ID 
1M15) [48]. Carbon = green, nitrogen = dark blue, oxygen = red and 
phosphorus = orange. In AK, Arg 28 o contacts the 0 3p oxygen of ADP, an 
oc-oxygen and the Asp 32 4 side chain. In ATP, the 0 3p oxygen bridges to 
the y-phosphate which is mimicked by nitrate in this transition state 
analog complex. Hydrogen bonds are shown with red dotted lines. 
doi:1 0.1 371 /journal. pone.01 0831 0.gOl 2 

G-proteins and other enzymes that have been discussed above. 
Thus, one should consider whether there could be a holistic 
rationalization of 0 3 p interactions, common to phosphagen 
kinases, G proteins and other enzymes. 

The 100-fold rate reduction of R280K is commensurate with 
changes on conservative mutations of arginine kinase Glu 2 25 and 
Gys27i that are implicated respectively in base catalysis and 
electrostatic modulation of the guanidino substrate [45,47]. Thus, 
the R280K mutant adds to the picture of phosphagen kinases 
employing several means of enhancing rate, of which the 0 3 p 
interactions are an important, but not dominant component. 
Furthermore, it is noted that the wild-type turnover is not limited 
by chemistry, but by protein conformational change [17]. Thus 
any of these mutations could be having a greater effect upon the 
chemistry, because we are only able to measure the part of slowing 
beyond the 104 s 1 WT turnover. 

Conclusions 

Several messages can be taken from the current work. Firstly, 
top-down analysis of protein structure databases can reveal 
commonalities in the disposition of diverse active sites catalyzing 
analogous reactions. Secondly, 0 3 p interactions figure prominent- 
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